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Abstract 


An unstaggered cascade consisting of only three syinnietrical airfoils (NACA 0015) was 
fabricated for experimental purpose. Wind tunnel tests were carried out to examine the 
aerodynamic characteristics of the airfoils in this configuration. Pressure distribution over 
the cascade airfoils were measured for different angles of attack with different space chord 
ratios and for two different Reynolds Numbers(Re.Nos : 3.31 x 10^, 4.6 x- lO*’). With 
the help of pressure plots, lift and drag were calculated. 


Wake measurements were made by means of wake survey rake a,ri(l total drag was calcu- 
lated. From the experimental studies, it was found that due to cascade effect the stall was 
delayed and some lift was obtained upto angle of attack 25°. 


The Cl vs a curve for an airfoil in the middle of the cascade has lower slope, compared 
to a single airfoil. As the angle of attack increases drag also increases, upto angle of attack 
15°. After that the drag decreases due to cascade effect. Thus the cascade effect improves 
lift-drag ratio for angles of attack above 15°. 
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Material Constants 
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Flow quantities 
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coefficient of viscosity of air 
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static pressure at the wake region 

total pressure at the wake region 

dynamic pressure 
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out flow velocity 

velocity at wake region 

Reynolds Number 

angular velocity of the turbine 


space chord ratio 
chord length 
wing area 
wing span 

gap between planes of the biplane turbine 
angle between planes of the biplane turbine 
l)itch between two blades(see Fig. 1.3) 
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Chapter 1 

INTRODUCTION 


1.1 Turbomachinery - Problem Outlook 

Two dimensional cascade tests and analysis have historically played a key role in the devel- 
opment of axial-flow compressor , axial-flow turbine and pump design. Various organizations 
contributed numerous cascade tests from which a wealth of design informations can be ob- 
tained. A successful investigation of the cascade problem needs systematic experimental 
research, as it was done for the single airfoil. In this work we wish to investigate the aero- 
dynamic performance of the symmetrical airfoil in unstaggered cascade formation. The idea 
of the present research work is that a deeper knowledge of the complex flow phenomena, as 
they occur in turbomachines cannot be acquired without experimental research. The first 
step in the investigation of the flow through cascades is the incompressible flow through a 
two dimensional cascade. The main problem in investigation of the two dimensional cascade 
is to find the relation between the geometric and aerodynamic parameters of the cascade. 
The geometric parameters are the profile of the Irlade with the maximum camber, leading 
edge radius, thickness ratio etc. The geometry of the cascade is given by space chord ratio, 
angle of stagger P etc. Present experiment d(!als with symmetrical airfoils, with thickness 
ratio t/c = 15 %, in unstaggered cascade formation with dilfcrent space chord ratios. These 
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parameters define the geometry of the cixscach; depe.iuling on the; arrangcmients of the blades 
in the cascade. 

The aerodynamic parameters of the cascade are 

inflow and outflow velocity, 

and angle of inflow and angle of outflow , 

deflection of the flow deviation , 

pressure difference across the cascade , 

chord wise pressure distribution over the blade, 

resultant aerodynamic force on the blade 

the tangential component of force , and the axial component of force, 
loss of total pressure in the cascade. 

The following parameters can be set at any required value within appropriate range. 

• inflow angle 

• inflow velocity 

The following quantities are measured as the function of these parameters, 

• 1. pressure distribution on the blade 

• 2. the total pressure behind the cascade 


1.2 Application To Wave Energy — A brief Outlook 

One of the most promising way of utilizing sea wave is provided by the oscillating water col- 
umn(OWC) wave energy device which is pollutant free and economically viable. The OWC 
device provides the simplest and possibly the most reliable means of converting irregular 
wave motion into high speed rotational movement required for electric power generation. For 
countries surrounded by the sea such as the British Isles, Japan, India, Portugal, Australia 
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and other island countries, wave energy conversion is an attractive technological proposition. 
The United Kingdom and Japan are pioneers in developing the wave energy and leading sev- 
eral activities. 


1.3 The Wells Turbine 

In the oscillating water column(OWC) device, the flow of air, displaced by the free surface 
within a chamber open at a immersed bottom, drives a turbine. The requirement of trans- 
forming the reciprocating airflow into unidirectional rotational motion is met without the 
need for any rectifying valves, by a turbine invented by Wells(1976). 

The Wells turbine is a self rectifying axial flow air turbine suitable for extracting the 
energy from reversed cyclic airflows. These conditions are encountered in the OWC caisson 
of wave energy converters. The Wells air turbine consists of several symmetrical airfoil blades 
set around a central hub as shown in Fig. 1.1 with their chord planes oriented normal to the 
axis of rotation. 

1.3.1 Principle of operation 

According to classical airfoil theory, an airfoil, which is s(J, at an angle of incidence alpha 
in a fluid flow, generates a lift force L, normal to the free stream and also a drag force D 
which is in the direction of the free stream. These lift and drag force can be resolved into 
tangential and axial components Fr ,Fx respectively. The. tangemtial and axial components 
are expressed as follows and also shown in Fig. 1.2 
Ft = L sin a — D cos a 


Fx = L cos a + D sin a 
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The tangential and axial force coefficients are Cpr Cp^ for a 2D airfoil. The force 
Ft is responsible for the power output, the force Fx results in an axial thrust along the axis 
of the rotor which has to be absorbed by the bearing. The turbine can start from rest (a 
= 90 degree) and accelerate to a running condition if the solidity of blades is higher than a 
certain minimum value. 

For symmetrical airfoil, the direction of tangential force Ft is the same for the both 
positive and negative values of 'a' . If such airfoil blades are positioned around an axis of 
rotation, this will always rotate in a single direction regardless of the direction of airflow. 
Therefore the turbine does not need rectifying valves in a bi-directional flow, the Wells 
turbine itself is a self rectifying turbine. That is why Wells wind turbine is well suited for 
wave energy conversion from OWC devices. 

There are two types of Wells turbine, mono-plane Wells turbine and bi-plane Wells tur- 
bine. Mono-plane turbine consists of symmetrical airfoil blades set around a central hub and 
bi-plane Wells turbine consists of two such rows of blades separated by a gap(Fig 1.3). One 
mono-plane Wells turbine is installed at Trivandrum in India. (Raghunathan 1995) 

1.4 The Cascade Wells Turbine 

Ghosh(1996a) has suggested a cascade turbine for the OWC wave energy device. This 
turbine is the subject of the present study.Tlie two view of the i)resent cascade Wells turbine 
and Wells turbine are shown in Fig. 1.4. Cascade turbine consists of two cascades at the 
top and bottom of the nacelle(Fig.l.4b) Each cascade consists of 6 blades lined up in the 
axial direction as in Fig.l.4b Whereas a typical Wells turbine consists of 12 blades in a single 
circular row. Chosh(1996a) has shown that the Working principle is same as VAWT(Vertical 
Axis Wind Turbine). Moreover the cascade effect allows a large angle of attack without stall 
in five upstream rows. 
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1.5 Objective of the Present Work 

From the past studies on Gas turbines ( compressor or turbine cascade) Pressure distribution 
of cambered airfoil cascades are available. From studies on Wells turbine, pressure distri- 
bution of tandem cascades of symmetrical airfoils is also available. Tandem cascade is also 
known as a cascade of 90 degree stagger angle. But little information about cascade of zero 
stagger symmetrical airfoils is available . So Present work is related to cascade studies of 
symmetrical airfoils. Schematic diagram of present cascade are shown in Fig 1.4. 

The aerodynamic parameters of the cascade are 
inflow and outflow velocity V'oo , Vi, 
and angle of attack o:, 

S/C space chord ratio, 

pressure difference across the cascade 6p — — Poo ^ 

chord wise pressure distribution over the blade p(x), 
resultant aerodynamic force on the blade L,D, 
loss of total pressure in the cascade Spo, 

The following parameters can be set at any required value within appropriate range. 

• inflow angle = a 

• inflow velocity = 1^0 

The aim is 

• to find the pressure distribution around the cascade airfoils with different space chord 
ratios at different angles of attack and Reynolds Numbers, 

• to find optimum space chord ratio needed for achieving high L:D, for angles of attack 
from 15° to 27.5°, 


• to find the total pressure behind the cascade 
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• to find the total drag which is calculated by wake survey method, 

• to find the separation points at different angles of attack and different space chord 


ratios for cascade airfoils. 



Chapter 2 

LITERATURE SURVEY 


The theory of the flow past cascade or infinite series of airfoils, regularly spaced is of im- 
portance in connection with the design of compressor, turbine and pump blades. Large 
amount of low speed cascade and high speed cascade data are available about cambered 
airfoil. Very little amount of experimental data are available about unstaggered cascade of 
symmetric airfoils. Schlichting(1954) made some investigations on unstaggered symmetric 
airfoil cascade. 

One of the most popular device in wave energy conversion is the water to air junction 
which is called oscillating water columns proposed by Masuda in .Japan. The popularity 
of the oscillating water column is derived from this simplicity which is a good feature of 
engineering design(Salter 1989). Dr. A. A. Wells invented a new design of turbine (which is 
called Wells turbine) at the Queen’s University of Belfest. A practical problem encountered 
with the Wells turbine is its difficulty to run up to operational speed when started from 
rest. This phenomenon where the turbine accelerates up to a certain speed which is much 
lower than the operational speed, is termed crawling. In order to avoid crawling , a high 
solidity and a low hub to tip ratio rotor is therefore required (Raghunathan 1982). Typically 
a hub to tip ratio of 0.6 and a solidity of 0.6 is required to overcome crawling. The blades 
in tandem cascade have better starting characteristics than an isolated blade, but will have 
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a lower ultiiiial.e cllicieney a(, I.1 h! oi)eial,ing c{)n(li(,ioii(ll.!i,gliuaa(,han 1981). 

There are many reports which describe the performance of tlie Wells turbine both at 
starting and at running conditions. According to their results, the wells turbine has inherent 
disadvantage, that is, lower efficiency, poor starting and a high axial thrust in comparison 
with conventional turbines. In order to overcome tlu^se weak points, an imi)ulse turbine with 
self-pitch-controlled guide vanes for wave energy conversion has been proposed(Kim et al 
1988). 

Ghosh(1996a , 1996b) suggests two novel cascade configurations of the wind turbine of the 
oscillating water column device. A probable configuration for guide vanes has been discussed 
(Ghoshl996b) by him. He suggested an alternative novel configuration which has a stator 
cascade in the middle with rotor cascades on either side. This turbine is unlikely to stall even 
in high velocity winds created by large waves. It operates at a relatively low tip speed ratio. 
However, the object of the study here is not this cascade, but a symmetrical airfoil cascade 
for a turbine of high tip speed ratio. (Ghoshl996a) It is argued that the cascade effect would 
delay onset of stall in the five upstream blades. 



Chapter 3 

EXPERIMENTAL SETUP AND 
PROCEDURE 


For the purpose of present study, we selected NACA 0015 syinmetrical airfoil. For this 
testing configuration, the dimension of the blades were selected according to the available 
space of the wind tunnel test section. The cascade is tested in closed circuit 2D wind tunnel. 

3.1 Fabrication of Airfoils 

Present study deals with cascade aerodynamics, so more than two airfoils were needed to 
do the experiment. Three airfoil blades were fabricated in Aero Workshop, IIT, Kanpur. 
The template was made, and used for checking the airfoil profile(NACA 0015). All the 
three blades were made out of wood. The span-wise location of pressure holes station were 
same for the all three blades and the station distance is 2Mnches from one end of the 
airfoil. Stainless steel tubes having 1.0 mm outer diameter and Ivini inner diameter were 
selected for pressure measurements. The bottom and top airfoils have same number of 
pressure tubes.(Fig.3.1)&(Tabular. 3.1) Middle airfoil has more number of pressure tubes 
because this airfoil has the recpiired cascade effect. Tlniy w('.r(^ i)laced at different chord- 


9 



10 


wise distances on both upper and lower surfaces of the, airfoil. A lar{>,(^ number of ijressure 
tubes were placed on the upper surface and fewer pressure tubes were placed on the lower 
surface. The distribution of pressure tubes are at close intervals near leading edge of the 
airfoil because pressure variation is rapid at this place. The remaining portion of the airfoil 
which has a small curvature, has a lesser number of jjressure tubes. Each pressure tube has 
length of Ziiiches and was fixed onto the groove by araldite. It should be noted that the 
airfoil blades were 2 dimensional, so there was no iumhI to fix longer length of presure tube 
on to the groove of the airfoil. After that, pressure holes were made on all tubes by Irm??, 
diameter drill bit. Finally, the surface of airfoil w(’re smootlnunul and i)olish(Hl by varnish. 
Details of dimension of blades and pressure tubes locations are shown in Fig 3.1 

3.2 Circular Model Base Disc 

The circular disc, made of wood, has dimensions of 20 inch dia and 1 inch thickness. A 
rectangular piece of wood linches x l^inches" was cut out from the center of the disc. In 
its place two MS metal strips were fixed, which were used for changing the space between 
the two airfoils. The circular disc was used for mounting the blades to simulate the hub 
effect on the test model. The model with base shown in Fig 3.2. (photograph). A schematic 
diagram of Expeimental Set up is shown in Fig.3.3. 

3.3 Wind Tunnel 

The experimental test were c.onducted in 2D arm of a 5D closc'.d-circuit Wind Tunnel in Low 
Speed Aerodynamics Laboratory at IIT, Kanpur. The low speed tunnel is a large venturi 
where the airflow is driven by two centrally located fans each connected to 15 hp motor 
drive. The schematic diagram of the 5D wind tumu'.l is shown in Fig 3.4. The airflow with 
pressure pi enters into the nozzle at velocity Vi, and area is Ax. The nozzle converges to 
a smaller area at the test section, where the velocity increases to Vco^ and the pressure 
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decreases to poo- After flowing over the model, the air is passed into a diverging duct (diffuser), 
where the area increases to A^, the velocity decreases to V3 and the pressure increases pa . 
Since , it is a closed circuit tunnel, the air from the exhaust is returned directly to the front 
of the tunnel. The tunnel test section dimensions are given below. 

• Test section Width = 1 ft 

• Test section Height = 4 ft 

• Test section Length = 5| ft 

3.4 Wake Survey Rake 

Pressure distribution is used to determine the lift , ])r(!s'sure drag and moment coefficients 
for the airfoils that are tested in a wind tunnel. Pressure distribution does not give the 
viscous drag. The total drag coellicient , in contrast, is usually obtained by measuring the 
momentum loss(or energy) of the air behind the airfoil. Data for this are obtained by a 
pressure-measuring device called a wake survey rake. Wake rake has a profile of symmetrical 
airfoil and made out of wood. It has 35 pitot probes and 5 static probes fixed equally spaced. 
Diagram shown in Fig 3.5(photograph) and Fig. 3. 3. Tlu; distance b(!twe<ui trailing edge of 
the airfoil and total pressure probe of wake rake is 4.3 C. 

3.5 Digital Micromanometer 

In this experiment Furness controls limited FC012 pressure scanner with two Digital Micro- 
manometers were used for measuring the model static pressure. It has a selection box with 
60 channels. A digital micromanomcter has a maximum range of ±200 imn of H-zO. 

The stainless steel tubes of the blade model were connected to FC012 pressure scanner 
through plastic tubes. The static pressure at every location was noted down just by selecting 
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:he channel in pressure scanner. This scanned ])r(’.ssur('. value sliows in digital rnicrornanorne- 
i;er. The micromanometer gives the difference between model surface static pressure and 
Free stream static pressure in iiviih of H-iO. I’lu^ pil.ot-sta.tic: probe was c;onnec;ted to another 
digital micromanometer which displayed directly the speed in m/s or difference between to- 
tal pressure and free stream static pressure in mm of H 2 O. Digital Micromanometer with 
Pressure Scanner are shown in Fig.3. 6. (photograph) and Fig.3.3. 

3.6 Water Manometer 

Water multimanometer has 40 tubes which were used for measuring the wake total and static 
pressure. 5 tubes were connected to static pre.ssure tubes and 35 tubes were connected to 
total pressure tubes of wake rake. Water manometer givc;s the pressure values in ynm. 

Before conducting the experiment the pressure holes- were checked for any blockage by 
passing compressed air. In the first set of readings, we found that two pressure tubes were 
blocked in the pressure scanner and these were cleared. The Wind tunnel was run at differ- 
ent speeds, as decided and the measurements were done for the different angles of attack of 
interest and different space chord ratios. Water multimanometer is shown in Fig 3.7(photo- 
graph) . 


3.7 Experimental Procedure 

3.7.1 Calculation of coefficient of Pressure 

Coefficient of pressure is defined as 

^ Psi. ~ Poo Psi Poo 

0.5 p 1 Poco Poo 

Where 

Psi - the static pressure of airfoil at particular pressure tube, 
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Poo " frse stream static pressurci, 

Pooo - free stream total pressure, 

In these experiments , static lu!a(l(measurecl with rcihucnc.e to ) and dynamic head 
were measured separately. So the ratio between static and dynamic head is coefficient of 
pressure. Coefficient of pressure are plotted vs X/C. We found that we did not get the 
reading for exact stagnation point (where Cp is +1) in many cases. To overcome this, the 
location of stagnation point in lower surface of the airfoil was taken from Shevell(1989) for 
NACA 0012. upto a — 10". For greater values of n, this stagnation point is located by 
extrapolating the curve. This is the defect of this experiment. Hence larger number of 
pressure tubes are needed in the lower surface of tln^ airfoil, for further experiments. 

3.7.2 Calculation of Lift and Drag from Pressure Measurement 

The aerodynamic forces on the airfoils are due to only two basic sources( Anderson 1991), 

1. Pressure distribution over the airfoil p(x), 

2. Shear stress distribution over the airfoil r 

The net effect of p(x) and r distril)ution inl.egral.cHl over the complete body surface is a 
resultant of aerodynamic force R and it can be split up into two components (sets of which 
are shown in Fig 3.8.) L is the lift component of R perpendicular to Poo, D is the drag 
component of R parallel to I4o , some times R is split into components perpendicular & 
parallel to chord C. 

F)v ■ normal force component of R perpendicular to chord C, 

Fc - axial force component of R paralhd to chord Cl, 

The geometrical relation between L, and D, Fc ar(' 

L = F/v t!ns a — Fc sin cv 


D = Fj\! sin fv + Fc cos n- 


(3.1) 

(3.2) 
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We note that 0 in Fig 3.9 i.s givc'.n by Uui ' {dy/dx) ainl .w is the. distane.e. along the airfoil 
curvature. The total normal and axial forces per , unit si)an are obtained 


Fn — - I iPu cos 0 + Tu sill 0) dsu 4- / [jH cos 0 - ri sin 0) dsi 
Jo Jo 


(3.3) 


rC rC 

Fc= {Pu sill 0 + Tu sin 0) ds^ + — {pi siu 6 + n cos 0) dsi (3.4) 

Jo Jo 

In these experiments, we measured static iiressure only, the shear stress term will be 
zero because we have neither measured nor calculated shear stress r The equation 3.3 & 3.4 
becomes. 


F!, = - 

- 1 {pu cos 0) dsu + , 

/ (/liCos 0)dsi 

(3.6) 


Jo j 

'o 



rC ri 

r; 

(3.6) 

Fc^ 

/ {pu sin 0) dsu - / 
Jo Jo 

(pfsin 0)dfii 


Eqn 3.5 & 3.6 can be written in nondimensional form. From the geometry shown in Fig 


0 = tan ^ {dy/dx) Ax = As x cos 0 

(3.7) 


Ay = {As X sin 0) 

(3.8) 


Area S = C x {1) 

(3.9) 

substituting eqn 3.7 & 3.8 into 3.5 & 3.6 and dividing by r/oo S and using eqn 3.9. We 

obtain the following integral forms for force coellicieiit 



X 

- 

X 

t> 

1 1 

(3.10) 

r - ^ 

i:(c;„% xAr,-tA'^ 

(3.11) 





The lift and drag coellicieiit (^iii b(! ohtaiiKul from ('(imition 3.1 & 3.2 


Cl = Ci.-fj COR a—C'i.-f. sin « 
Cd — Cpfj sin a + cos o; 

Cl is the lift coefficient 

Cd is the pressure drag coefficient 


(3.12) 


1.8 Calculation of Drag by the Wake Survey Method 

'he total drag of the cascade drag may be obtained by comparing the momentum in the air 
head of the model with momentum behind the model (Alen Pope 1947). 


D = X change in velocity 

sec 

D = 1 1 p (v„ K, - Cj) da 


Where, 

Voo is the free stream velocity, 

V'oo)is the wake region velocity, 

da is the small area of tlu'. wa.k(! iierjxuidicular to sl.itiain, 


(3.13) 

(3.14) 


Also, 



D 

0.5 p S 


Cd = 2 



[-2 


1/2 
* oo 


da 

'S 


(3.16) 


V'c 


2 Qoo 


y 

» Hi 


2 (ho 


P 


P 


(3.16) 



For a unit section of the airfoil S = C xl, and. the area da — dy x 1 , where y is measured 
■pendicular to the free stream. 



doo j ^ 


( 3 . 18 ) 


The total drag is C,i, = Cdp + 6^,;., . Where, 

Cdp is the pressure drag, is tlie skin friction drag. 



Chapter 4 


RESULTS AND DISCUSSIONS 

4.1 Cp Distribution of Cascade Airfoil 

The surface pressure measureiiieiit for tlie cascade airfoils were done by means of digital 
micromanometer for angles of attack of 0" to 27.5", and S/C = 0.5, 0.7, 0.9 for two different 
Reynolds Numbers(i?e No = 3.31 x 10^ at Coo = 25 in/s, Re No = 4.6 x 10® at Coo 
= 35 m/s). 

4.1.1 Cp for S/C = 0.5, a = 0^ to 27.5", at Uoo = 25 m/s & 35 m/s. 

The Cp distribution curve could not be closed in most cases and area under the curve was 
compared by eye estimation. 

• In Cp .distribution, shown in Figs. 4.1 & 4.2 it is seen that top airfoil upper surface 
nearly behaves like a single airfoil of NACA ()015(Abbott et al(1959)). Bottom airfoil 
lower surface shows slightly different pressure from a single airfoil. There may be half 
a degree error in setting of this airfoil. Tlui Cp of the lower surface of top airfoil and 
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ths upper surface of bottom airfoil are interfered with due to cascade effect and they 
have greater acceleration of the flow, manifested in greater negative pressures. The 
two surfaces of middle airfoil have almost identical pressure distribution, that means 
its producing almost zero lift at an angle of attack 0" . Also we perceive greater 
acceleration of the flow on both the surfaces compared to a single airfoil. This is due 
to cascade effect. 

• Figs 4.3 & 4.4, the separation occurs on upper surface of the top airfoil near trailing 
edge(TE). Bottom airfoil produces more positive lift compared to other airfoils in 
cascade. This is to be expected since separation is suppressed on the upper surface 
of the bottom airfoil due to cascade effect. The upper surface of middle & bottom 
airfoils are comparable. But the lower of the middle has more accelerated flow due to 
cascade effect. Hence middle airfoil produces less lift than bottom airfoil. For the sake 
of convenienc.e from now on we shall refer to the airfoils as top af, middle af & bottom 
af. 

• The top af in Figs. 4.3 &:4.4 produces almost no lift, since the lower surface has more 
than usual acceleration of flow due to cascade (dfect. 

• In Figs 4.3 to 4.10, the separation point gradually moves towards leading edge(LE) 
of top af, as the angle of attack increases from 0" to 17.5". And lower surface of 
the middle af have positive value of Cp. The Cp of lower surface of bottom af are 
consistently larger comi)ared to other af lower surfaces. Bc'.twe.en n'15"&17.5", the lift 
(area between upper & lower surface curves) of top af seems to reach a maxima. This 
signifies onset of stall. 
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• At high angles of attack, 20" in Fig 4.11(a), top af has stalled completely at Voo= 25 
m/s . Considering upper surface alone in Fig 4.11(a) , it can be seen that pressure 
recovery has almost stopped at the leading edge. But at = 35 m/s in Fig 4.12(a) 
the same upper surface. has slight pressure recovery. This may be due to increased 
Reynolds Number. The same trend is happening upto a — 22.5" . But at high angles 
of attack, Cp of upper surfac.es of middle af and bottom af have higher negative values. 
This is because of cascade effect. 

4.1.2 Cp for S/C = 0.7, cr = 15" to 27.5", at Foo = 25 m/s & 35 m/s. 

• Considering Fig 4.19(b), upper surface of middle af has negative pre,ssure of greater 
magnitude when compared to Fig 4.7(b). This is because of increase in S/C (S/C 
= 0.7). That means that the How accelerat(!S more , and middle af i)roduces more 
lift. Comparing Fig 4.19(a) and 4.7(a), the lift of both are same. The separation 
point occurs near the mid chord. In Fig 4.20, the upper surface of the top af has 
same pressure recovery as Fig 4.19. The upper surface of middle af has lower negative 
pressure when compared to 4.19(b). Fig 4.21(a) top af at a — 17.5", at Voo — 25 m/s, 
upper surface Cp distribution is almost flat along the chord length. It seems that top 
af has stalled. Fig 4.21 bottom af produces more lift compared to other af in cascade. 
In Fig 4.22(a) V"oo = 35 ni/s upp(!r surface! n(’ga,tivc pressure! is high. 'I'lds may be due 
to increased Reynolds Number. Fig 4.22(c) bottom af upper surface negative pressure 
is high at LE when compared te) Fig 4.21(c.). This may be due to increased Reynolds 
Number. Fig 4.23(a) top af has stalled completely at 14o = 25 ni/s. But in Fig 4.24(a) 
top af , upper surface negative pressure is high in a small zone near the leading edge. 
This is due to increased Reynolds Numb(!r. In Fig 4.24(b) middle af ui)per surface has 
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lower negative pressure compared to Fig 4.23(b). This is due to lesser flow acceleration. 

• At high angles of attack, cx = 20" to 27.5", the same trend happens. One thing should 
be noted that negative pressure magnitude of the upper surface of the middle af is 
increased near leading edge, when compared to S/C = 0.5, ct = 17.5" to 27.5" (Fig 
4.21 to Fig 4.30). This is because of increasing S/C ratio and it behaves a little like a 
single af. 

4.1.3 Cp for S/C = 0.9, a — 15° to 27.5°, at Voo — 25 m/s Sz 35 m/s. 

• For angles of attack from 15" to 20", similar trends as for S/C - 0.7 are observed. 

• Fig 4.35 a = 20" , top af has stalled and middle and bottom af pressure distribution 
has higher negative on the upper surface. But in Fig 4.36(a) upper surface of top af 
pressure is starting to decrease and then recover, this indicates that stalling is not total 
yet. The same trend is happening upto o; = 22.5°. When a is increased to 25" top 
af(Fig 4.39) stalls completely. The middle and bottom af also begin to stall. In Fig 
4.39, the middle af and bottom af have also stalled in the sense that the upper surface 
flow separates at the leading edge. But even then the cascade effect brings some benefit 
and keeps the upper surface pressure coefficient at a higher negative, when compared 
to that of top af. Thus the middle and bottom af j)roduce some! lift and less drag. The 
same trend continues from Fig 4.39 to Fig 4.42. 


4.2 Velocity Profile from Wake Measurement 

Wake region velocity profile vs location of total pressure prolies in the wake rake are plotted 
for all the angles of attack(cv = 0" to 27.5" ). Bezier approximation curve was used to get 
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smooth velocity profile. These experiments were done in 2D closed circuit wind tunnel, 
ike rake was mounted in test section which is used for single af wake measurement, but 
isent experiment is about cascade afs. Wake rake height is not quite sufficient for these 
periments. 

2.1 Velocity profile for S/C = 0.5, a = (D to TI Xf 

• In Fig 4.43 (a) & (b), at cv = 0° , the maximum loss of velocity is at region of top af 
and bottom af. When the angle of attack is increased upto 15” the loss of velocity is 
more behind the bottom af when compared to the other afs(Fig 4.44 to Fig 4.47). This 
is because more upstream flow seems to be deflected towards the top of the cascade. 
When a = 20” , the momentum loss is more behind top af. This is because of top af 
stalls. The momentum loss is lesser for middle af. This happens ui)to angles of attack 
22.5” (Figs 4.48 & 4.49). In Fig 4.50 to Fig 4.51 beyond rv = 25”, the trend is entirely 
different, that is middle af has a greater momentum defect. 

2.2 Velocity profile for S/C = 0.7, a = 15® to 27.5® 

• The velocity profiles are shown in Fig 4.52 to Fig 4.57. At Voo is 25 m/s the loss of 
velocity is more towards the top af than bottom af. In Fig 4.52(b) top and bottom 
af has almost equal amount of momentum loss. In fig 4.53(a), at Voo is 25 m/s the 
top af wake is lossing more momentum. Consider Fig 4.53(b) the maximum loss of 
momentum is in between toi) and middle afs at cv — 17.5". Fig 4.54(a) is not correct, 
may be this is due to experimental error. The momentum defect is increased at high 
angles of attack like a = 20” to 27.5" ( Fig 4.55 to Fig 4.57). When compared to 
Fig 4.43 to Fig 4.52, it seems that velocity loss in the case of S/C — 0.7 is more. This 
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because of increased blockage of the tunnel due to increased width of the cascade. 

2.3 Velocity Profile for S/C = 0.9, a = 15^^ to 27.5" 

• Considering Fig 4.58 (a) , the momentum defect has a single maximum in wake region. 
Fig 4.58(b) is not correct , may be this is because of experimental error. In Fig 4.59 
Fig 4.60, Fig 4.61, momentum (hdect is larg(i. Velocity prolih; se(uns to be a i)arabolic 
curve. At high angles of attack like a — 25" & 27.5" the momentum defect is more 
between top af and middle af. 

2.4 Drag Calculated by Wake Survey Method 

• Drag calculated by wake survey method , in Figs 4.64 and 4.65 are plotted from the 

experimental result. Drag increases when the angle of attack increases from 0^ to 15^, 

after that drag decreases when a increases for l oo = h I oo = 357u/s, 

S/C = 0.5. For NACA 0015 single af , drag increa.se.s with a (Critzes et al (1955)). 
When comparing single af and present cascade afs, the drag is lesser in cascade afs for 
higher angles of attack. We have measured wak('. measurement for S/C = 0.7 & 0.9 at 
, a = 15" to 27.5" with two different Reynolds Numbers. From Fig 4.64 and Fig 4.65 , 
we see that same trend is happening like S/C = 0.5. But with different drag values. 


.3 Cl VS a 

• Cl VQ a are plotted only for S/C = 0.5, a = 0,5, 10 degrees. For these angles of 
attack, Cp distribution has a closed curve(Fig 4.66 (a),(b),(c)). The Ci vs a are plotted 
in different graphs for top af , middle af and bottom af. Toj) af gives negative lift at 
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q; = 0°. After a is increased lift is also increasing . As the velocity is increased from 
25 m/s to 35 m/s, the lift also increases. Middle, af produce positive lift and bottom 
af is also producing more positive lift compared to othcu- afs. 

:.4 Conclusion 

he following conclusion were obtained, 

• In unstaggered symmetrical cascade af stall was delayed & some lift was obtained upto 
Q!= 25° 

• The Cl vs a curves for middle af & bottom af have lowcir slope, compared to a single 
af. 

• As the angle of attack increases drag also increases upto o; = 15°, after that the drag 
decreases. 


C.5 Further Studies 

• The same experiments should be performed with more luessure tubes near the LE of 
the af and lower surface of all the afs. So that Cj, distribution curve can be closed and 
one can calculate the lift and drag. 
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Figun; i.l: Sclicniulac of Wcdls. Tiirbiiu'. 








Figure 1.2: Tangential and Axial foree.s of Turbine Blade 





Figure 1.3: The biplane Wells Turbine 







Figure 3.2: The Model with Base Disc bid'ore inouiitiiig in tlie Test Section 
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Figure 3.5: Wake rake in the Wind Tunnel Test Section 



Figure 3.6: Digital Micromanoineter with Pressure Scanner 



Figure 3.7: Water Multiiriaiiuiueter 








x/c 

Tube No 

X/C 

Tube No 

0 

19 

0 

39 

0.00714 

20 

0.00714 

40 

.37 

21 

0.375 

41 

0.643 

22 

0.643 

42 

0.821 

23 

0.821 

43 

0.005.4 

24 

0.0054 

44 

0.107 

25 

0.107 

45 

0.161 

26 

0.161 

46 

0.214 

27 

0.214 

47 

0.259 

28 

0.259 

48 

0.375 

29 

0.304 

49 

0.5 

30 

0,375 

50 

0.571 

31 

0.446 

51 

0.643 

32 

0.5 

52 

0.714 

33 

0.574 

53 

0.786 

34 

0.643 

54 

0.857 

35 

0.714 

55 

0.929 

36 

0.786 

56 


37 

0.857 



38 

0.929 
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0.00714 

0.375 

0.643 

0.821 

0.0054 

0.107 

0.161 

0.214 

0.259 

0.375 

0.5 

0.571 

0.643 

0.714 

0.786 

0.857 

0.929 


Table 3.1: Location of- Pressure Tubes 




(c) Bottom Airfoil 


Figure 4.1: Pressure Distribution of Airfoil; S/C — 0.5, cr — Vqo ~ 
* Abbott et al (1959), + Lower Surface, x Upper Surface. 








(c) Bottom Airfoil 


Figure 4.2: Pressure Distribution of Airfoil; S/C = 0.5, o; = 0", Poo 
* Abbott et al (1959), + Lower Surface, x Upper Surface, 









(c) Bottom Airfoil 


Figure 4.3: Pressure Distribution of Airfoil; S/C — 0.5, 

+ Lower Surface', x Upper Surface 









(c) Bottom Airfoil 


Figure 4.4: Pressure Distribution of Airfoil; S/C — 0.5, ( 

+ Lower Surface, x Upper Surface 
















(c) Bottom Airfoil 


Figure 4.6: Pressure Distribution of Airioil, S/C O.o, c 

+ Lower Surface, x Upper Surface 







(c) Bottom Airfoil 


Figure 4.7: Pressure Di,st.ril)ut.iuu of Airfoil; S/C == 0.5, r 

+ Lower Surface, x Upper Surface 









(c) Bottom Airfoil 


Figure 4.8: Pressure Distrihutiou of Airfoil; S/C = 0.5, n 

+ Lower Surface, x Upper Svirface. 







(c) Bottom Airfoil 


Figure 4.9: Pressure Disl;ril)utiou of Airfoil; S/G = 0.5, n 

+ Lower Surface, x Ui)per Surface 

















Figure 4.11; Prei 


(c) Bottom Airfoil 


ire 


Distribuliou of Airfoil; S/C = 0.5, t 
+ Lower Surface, x Uiiiier Stirface, 






















(c) Bottom Airfoil 


Figure 4.14: Pressure Distribution of Airfoil; S/C = 0.5, c 

+ Lower Surface, x Upper Surface 
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(c) Bottom Airfoil 


Figure 4.15: Pressure Distribution of Airfoil; S/C = 0.5, ( 

+ Lower Surfat;e, x Uirper Surface, 







(c) Bottom Airfoil 


Figure 4.16: Pressure Distribution of Airfoil; S/C = 0.5, 

+ Lower Surface, xUititer Surface 









(c) Bottom Airfoil 


Figure 4.17: Pressure Distribution of Airfoil; S/C = 0.5, r 

+ Lower Surface, x Upper Surface 








(c) Bottom Airfoil 


Figure 4.18: Pressure Distribution of Airfoil; S/C = 0.5, o 

+ Lowcu' Surface, x Ui)i)er Surface 








(c) Bottom Airfoil 


Figure 4.19: Pressure Distribution of Airfoil; S/C = 0.7, 

+ Lower Surface, x Ui)i)er Surface 







(c) Bottom Airfoil 


Figure 4.20: Pressure Distribution of Airfoil; S/C = 0.7, < 

+ Lower Surface, x Upper Surface, 









(c) Bottom Airfoil 


Figure 4.21: Pressure Dist,ril)Ut,i()U of Airfoil; S/C = 0.7, < 

+ Lower Surface, x Upper Surface 
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(c) Bottom Airfoil 

Figure 4.22: PreHsuro Di.si.rihiit.ion of Airfoil; S/( ! = 0.7, r 

+ Lower Siirfaee, x Upper Surface 











(c) Bottom Airfoil 


Figure 4.24: Pressure Distriloutioii of Airfoil; S/C i 

+ Lower Surface, x Ui)i)er Surface 

















(c) Bottom Airfoil 


Figure 4.26; Pressure Distribution of Airfoil; S/C = 0.7, c 

+ L(rw(!r Surface, x Uirpc'r Surfaci’ 









(c) Bottom Airfoil 


Figure 4.27: Pressure Distribution of Airfoil; S/C — 0.7, < 

+ Lower Surface, x Upper Surface, 















(c) Bottom Airfoil 


4.29: Pressure Distribution of Airfoil; S/C = 0.7, r 
4- lj()W('r Surface, x Upjx’r Surface 










(c) Bottom Airfoil 


Figure 4.30: Pressure Distribution of Airfoil; S/C = 0,7, o 

4- Lower Surface, x Ui)i)er Surfacie 
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(a) Top Airfoil 
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(b) Middle Airfoil 
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(c) Bottom Airfoil 


Figur6 4.31: PrGssurG Distribution of Airtoilj S/G — 0.9, cr 15 , Vqq — 25 ixi/s 

+ Lower Surface, x Ui)per Surface. 






(c) Bottom Airfoil 


Figure 4.32: Pressure Distribution of Airfoil; S/C = 0,9, 

+ Lower Surface, x Ui)per Surface 









(c) Bottom Airfoil 


Figure 4.33; Pre.ssure Distribution of Airfoil; S/C = 0.9, rv = 17.5", K 

+ Lower Surface, x Upper Surface. 








(c) Bottom Airfoil 


Figure 4.34; Pressure Distribution of Airfoil; S/C = 0.9, a = 17.5°, K 

+ Lower Surface, x Upper Surface. 








(c) Bottom Airfoil 


Figure 4.35: Pressure Distribution of Airfoil; S/C = 0.9, i 

+ Lower Surface, x Upper Surface 








(c) Bottom Airfoil 


Figure 4.36: Pressure Distribution of Airfoil; S/C — 0.9, a 

+ Lower Surface, x Uitpttr Surface. 









(c) Bottom Airfoil 


4.37: Pressure Distribution of Airfoil; S/C — 0.9, a — 22.5^^ 14 
+ Lower Surface, x U!)i)er Surface. 



























(c) Bottom Airfoil 


‘igure 4.40: Pressure Distribution of Airfoil; S/C 

+ Lower Surface, x Upper 









(c) Bottom Airfoil 


4.41: Pressure Distribution of Airfoil; S/C = 0.9, c 
+ Ijovvcr Surfaces, x Uirixu- Surface 
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Figure 4.43 Velocity Profile from Wake M('asuremeiit; S/C = 0,5, a = 0 ", 
oo = 25 m/s, (b). Voo = 35 m/s, + Experimental, - -- Bezier Apporoximation Curve, 
)cation of Top Airfoil, (DLocatioii of Middle Airfoil, 0 Location of Bottom Airfoil. 



Figure 4.44 Velocity Profile from Wake Measurement; S/C = 0.5, a = 5", 

Voo = 25 m/s, {b). V'oo = 35 m/s, 4 - Experimeiita.l, - -- Bezier Apporoximation Curve, 
Location of Top Airfoil, @ Location of Middle Airfoil, 0 Location of Bottom Airfoil. 
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Figure 4.4& Velocity Profile from Wake Meaaulemeut; S/C = 0.5, a = 10”, 

„ = 25 m/s, (b). Voo = 35 m/s, + Experimental, - -- Bezier Apporoximation Curve, 
cation of Top Airfoil, (g) Location of Middle Airfoil, @ Location of Bottom Airfoil. 
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Figure 4. Velocity Profile from Wake Mea,surement; S/C = 0.5, a = 15”, 

Voo = 25 m/s, (b). Voo = 35 m/s, + Experimental, - -- Bezier Apporoximation Curve, 
vocation of Top Airfoil, (?) Location of Middle Airfoil, 0 Location of Bottom Airfoil. 








gure 4. A7 Velocity Profile from Wake Measurement; S/C = 0.5, rv = 17.5", 

= 25 to/s, (5). Voo = 35 m/s, + Experimental, - Bezier Apporoximation Curve, 
tion of Top Airfoil, © Location of Middle Airfoil, (^Location of Bottom Airfoil. 



''igure 4. 4^ Velocity Profile from Wake Measurement; S/C = 0.5, « = 20", 

= 25 m/s, {b). Voo — 35 m/s, + Experimental, - -- Bezier Apporoximation Curve, 
ition of Top Airfoil, @ Location of Middle Airfoil, ^Location of Bottom Airfoil. 
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'igure Velocity Profile from Wake Measurement; S/C = 0.5, cv = 22.5", 

, = 25 m/s, {b). Voo = 35 m/s, + Experimental, -- - Bezier Apporoximation Curve, 
ation of Top Airfoil, ^ Location of Middle Airfoil, @ Location of Bottom Airfoil. 



Figure 4. Velocity Profile from Wake Measurement; S/C = 0.5, a = 25", 

'oo = 25 m/s, {b). Voo — 35 m/s, + Experimental, - -- Bezier Apporoximation Curve, 
Dcation of Top Airfoil, ^ Location of Middle Airfoil, 0 Location of Bottom Airfoil. 
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Figure 4*51 Velocity Prolile from Wake Measurement; S/C = 0.5, a = 27.5°, 
o = 25 ?u/s, (6). Koo = 35 m/s, + Experimental, - -- Bezier Apporoximation Curve, 
cation of Top Airfoil, © Lo(uition of Middle Airfoil, < 3 ) Location of Bottom Airfoil. 





'’igure 4.52; Velocity Profile from Wfikc' Measuif'ment; S/C = 0.7, (v = 15”, 

= 25 m/s, {b). V'oo = 35 rn/s, + Experiiiiental, - -- Bezier Apporoximation Curve, 
ation of Top Airfoil, ^Location of Middle Airfoil, & Location of Bottom Airfoil. 



Figure 4.53: Velocity Profile from Wake Measurement; S/C = 0.7, a — 17.5”, 

/x, = 25 m/s, (6). Voo = 35 m/s, + Experimental, - - - Bezier Apporoximation Curve, 
ocation of Top Airfoil, Location of Middle Airfoil, 0 Location of Bottom Airfoil. 







figure 4.54: Velocity Prolile from Wake Measurement; S/C = 0.7, a = 20", 

= 25 m/s, (6). Voo = 35 m/s, + Experimental, - Bezier Apporoxiination Curve, 
ation of Top Airfoil, 0 Location of Middle Airfoil, 0 Location of Bottom Airfoil. 



Figure 4.55: Velocity Profile from Wake Measurement; S/C = 0.7, a = 22.5", 

X, = 25 m/s, (b). Voo = 35 m/s, + Experimental, - -- Bezier Apporoxiination Curve, 
cation of Top Airfoil, 0 Location of Middle Airfoil, 0 Location of Bottom Airfoil. 
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’igure 4.56: Velocity Profile from Wake Measurement; S/C = 0.7, « = 25", 

= 25 m/s, {b). Voo = 35 m/s, + Experimental, Bezier Apporoximatioii Curve, 
ition of Top Airfoil, ® Location of Middle Airfoil, Location of Bottom Airfoil. 
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Figure 4.57: Velocity Profile from Wake Measurement; S/C = 0.7, a ~ 27.5", 

'oo = 25 m/s, (6). Voo = 35 m/s, + Experimental, - -- Bezier Apporoximatioii Curve, 
ication of Top Airfoil, (l)Location of Middle Airfoil, (^Location of Bottom Airfoil. 







gure 4.58: Velocity Profile from Wake Measurement; S/C = 0.9, a = 15°, 

= 25 m/s, {b). Voo = 35 m/s, 4- Experimental, - -- Bezier Apporoxiraation Curve, 
tion of Top Airfoil, Location of Middle Airfoil, (S) Location of Bottom Airfoil. 
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Figure 4.59: Velocity Profile from Wake Measurement; S/C = 0.9, a = 17.5°, 
oo = 25 m/s, {b). Voo = 35 rn/s, -f Experimental, - Bezier Apporoximation Curve, 
)cation of Top Airfoil, @ Location of Middle Airfoil, 0 Location of Bottom Airfoil. 






ire 4.60: Velocity Profile from Wake Measurement; S/C = 0.9, a = 20°, 

25 m/s, (b). Voo = 35 m/s, + Experimental, - -- Bezier Apporoximation Curve, 
m of Top Airfoil, Location of Middle Airfoil, Location of Bottom Airfoil. 




.figure 4.61: Velocity Profile from Wake Measurement; S/C = 0.9, a = 22.5°, 
o = 25 m/s, {b). Voo = 35 m/s, + Experimental, - -- Bezier Apporoximation Curve, 
cation of To]) Airfoil, (^Location of Middle Airfoil, Cl) Location of Bottom Airfoil. 





Figure 4.62: Velocity Profile from Wake Measiireinent; S/C = 0.9, a = 25®, 

, — 25 m/s, (b). Voo = 35 m/s, + Experimental, - Bezier Apporoxirnation Curve, 
:ation of Top Airfoil,^ Location of Middle Airfoil, ^Location of Bottom Airfoil. 
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Figure 4.63: Velocity Profile from Wake Measurenient; S/C = 0.9, a = 27.5®, 

Voo = 25 m/s, {b). Voo = 35 m/s, + Experimental, - -- Bezier Apporoxirnation Curve, 
Location of Top Airfoil, © Location of Middle Airfoil, 0 Location of Bottom Airfoil. 
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Figure 4.64: Drag vs a, \ = 25 in/; 
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pendix 


ier approximation curves 

curves and surfaces are credited to P. Bezier of the frencli car firm Regie Renault who 
ped it about (1962) and used them in his software system UNISERF, which has been 
»y designers to define the outer panels of several Renault cars. These curves are known 
der curves. 

athematically for n+1 control points, the Bezier .curve is defined for the following 
omial of degree n, 

■P(^) = ELo Pi Bi,n{u) 0< u< 1 (1) 

herer P(u) is any point on the curve and Pi is a control point. Bi,n are the Bernstein 
;s polynomials. Thus, the Bezier curve has a Bernstein basis. The Bernstein polynomial 
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blending or basis function for the Bezier curve and is given by 

Bi,n{u) = C{n,i) v} (1 - uY-^ (2) 

C(n,i) is the binomial coefficient 

C{n,i) = j (3) 

iig equation 2 and 3 and observing that C(n,o) = C(n,n) =1, equation 1 can be 
to give 

P{u) = Foil- n)” + Pi C{n, 1) «(1 - C{n, 2) u^(l - 

+ + Pn-i (l-u) + P„tP 0 < u< 1 (4) 

iharacteristics of Bezier curves are 

e basis functions are real. 

le degree of polynomial defining the curve segment is one less than the numbering 
defining polynomial points. 

hey generally follow the shape of the defining polygon. 

he first and last points of the curve are coincident with the first and last points of 
le defining polygon. 

'he curve is contained within the convex hull of the defining polygon. 

."he curve interpolates the first and last points. That is it passes through the first and 
ast points Po and Pn if we substitute u = 0 and 1 in eqn. 4 

The curve is tangent to the first and last segments of the characteristic polygon. 

3ference, 

Ibrahim Zeid(19991), CAD/CAM Theory and Practice., McGraw-Hill, Inc., New York. 



